On Indirect CP Violation and Implications for D° — D° and B s — B s mixing 
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The two kinds of indirect CP violation in neutral meson systems are related, in the absence of new 
weak phases in decay. The result is a model-independent expression relating CP violation in mixing, 
CP violation in the interference of decays with and without mixing, and the meson mass and width 
differences. It relates the semileptonic and time-dependent CP asymmetries; and CP-conjugate 
pairs of time-dependent D° CP asymmetries. CP violation in the interference of decays with and 
without mixing is related to the mixing parameters of relevance to model building: the off-diagonal 
mixing matrix elements |Pi2|, and <f>i2 = arg(Mi2/Pi2). Incorporating this relation into a fit 

to the D° — D° mixing data implies a level of sensitivity to \4 I ?2\ of 0.10 [rad] at la. The formalism 
is extended to include new weak phases in decay, and in Ti2 . The phases are highly constrained by 
direct CP violation measurements. Consequently, the bounds on \4>\2\ are not significantly altered, 
and the effects of new weak phases in decay could be difficult to observe at a high luminosity flavor 
factory (D°) or at the LHC (B s ) via violations of the above relations, unlike in direct CP violation. 
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I. INTRODUCTION 

There are two kinds of indirect CP violation in neutral 
meson decays, CP violation in pure mixing (CPVMIX) 
and CP violation in the interference of decays with and 
without mixing (CP VINT) (see, for example, [1, 2]). Let 
M° and M° be the interaction eigenstates of a neu- 
tral meson system. Indirect CP violation in pure mix- 
ing is due to a non-vanishing relative phase, </>i2 = 
arg(Mi 2 /ri 2 ), between the dispersive (Mi 2 ) and absorp- 
tive (ri2) parts of the M° — M° transition amplitude. It 
is responsible for CP asymmetries in semileptonic decays 
(M°, M° -> e ± X). Indirect CP violation in the interfer- 
ence of decays with and without mixing (M° — > M — » / 
and M° — > /) can occur in decays to final states which 
are common to M° and M°, leading to time-dependent 
CP asymmetries. 

Direct CP violation corresponds to different magni- 
tudes for decay amplitudes related by CP conjugation. It 
requires at least two amplitude contributions with differ- 
ent CP violating weak phases and different CP conserv- 
ing strong phases. The weak phases present in the decay 
amplitudes in addition to the dominant Standard Model 
(SM) weak phase, subsequently referred to as "new weak 
phases in decay" , can also lead to unequal CPVINT mea- 
surements for different final states, and to T-violating 
triple-product correlations for V V final states Q , even if 
strong phase differences are absent. 

In general, CPVINT receives contributions from 
CPVMIX and from new weak phases in decay. However, 
if the latter is absent, then CPVINT originates solely 
from the mixing phase </>i2, and therefore it must be con- 
nected to CPVMIX. Consequently, two related formulae 
can be derived: (i) an expression for CPVINT in terms of 
the mixing parameters </>i2, |Mia|, and |r i2 1 , see Eq. (f5*2"j). 
Such a relation was first derived in the limit Mi 2 <C T\2 
[H ; (ii) a model- independent expression relating the four 
mixing observables, i.e., the two kinds of indirect CP vio- 
lation and the neutral meson mass and width differences, 
see Eq. (IMl) . (i) allows a fit of the three mixing parame- 



ters to the four observables to be performed, (ii) leads to 
model-independent correlations between time-dependent 
and semileptonic CP asymmetries. It also leads to sim- 
ple relations between CP-conjugate time-dependent CP 
asymmetries in D° decays to non-CP eigenstates. 

Examples in which the connection between the indirect 
CP asymmetries can be realized are provided by the tree- 
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D° 
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level dominated decays, e.g., K° 
D° -> K + K~, tt+ti--, and B s -> J/*0. Contributions 
to these decays beyond the SM tree-level charged current 
interactions can be neglected in the SM itself, as well as 
in many of its proposed extensions. Thus, the underly- 
ing hypothesis of no new weak phases in decay is often 
valid. Of particular interest are applications to the D° 
and B s systems, where non- vanishing indirect CP asym- 
metries would constitute a clear signal for new physics. 
In many SM extensions they could be present at levels 
which can be measured at ongoing, imminent, or planned 
experiments. 

We review the neutral meson mixing and CP violation 
formalism in Sections II and III. Attention is paid to 
the independence of physical observables with respect to 
the sign convention for the neutral meson mass or width 
differences. In Section IV we derive the expression for 
CPVINT in terms of <j> 12 , \Mi 2 \, and |r 12 1 ; the model- 
independent relation between CPVMIX and CPVINT; 
and the resulting correlations among the time-dependent 
and semileptonic CP asymmetries. In the case of the D°, 
we discuss singly Cabibbo suppressed (SCS), Cabibbo 
favored (CF), and doubly Cabibbo suppressed decays to 
CP and non-CP eigenstates. In the case of the B s , we 
focus on b — > ccs mediated transitions, e.g., B s — > J/^(j>. 
In Section V a fit to the D° — D° mixing data is carried 
out to determine the allowed ranges for (f>f 2 , \M^\, and 
Irf^l in models with negligible new weak phases in the 
tree-level dominated decays. 

In Section VI we discuss in detail how the above re- 
sults would be modified by the appearance of subleading 
weak phases in the decay amplitudes, and in Ti2. Order 
of magnitude bounds on these weak phases, hence on vio- 
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lations of the relations between CPVMIX and CPVINT, 
can be obtained from existing direct CP violation mea- 
surements. It then follows that (i) the bounds on <ff 2 
(and |-Mj^|, |rf^|) do not change significantly, and (ii) it 
could be difficult to detect violations at currently allowed 
levels in the D° and B s systems, at a super B factory and 
at the LHC, respectively. In fact, the existence of new 
weak phases in decay would be much easier to discover 
directly, via direct CP asymmetry measurements. How- 
ever, with sufficient statistics it could be possible to iso- 
late and measure shifts in arg(r 12 ) (due to new physics, 
or to subleading 0(V* s V u b/V* s V c b) SM contributions in 
B a mixing) from such violations. We conclude in Section 
VII. 

While this work was in progress, Ref. [5[ appeared, 
which also explores the relation between the two indirect 
CP asymmetries, in the absence of new weak phases in 
decay. Our starting point for the derivation of a model- 
independent relation differs, in that it explicitly removes 
a discrete ambiguity in </> 12 <-* —4>i2, and allows us to ob- 
tain a simple general expression. The reader is referred 
to Q for a discussion of all four neutral meson systems. 
Also see @ for a discussion, based on 0], of correlations 
between time-dependent and semileptonic CP asymme- 
tries in decays to CP eigenstates. After completion of 
this work, we discovered that our model-independent re- 
lation, Eq. (|54[) . can be found in [7|. We augment their 
presentation by providing its derivation from the neutral 
meson mixing formalism, and by discussing its signifi- 
cance for relating CPVMIX and CPVINT. Our fit proce- 
dure for D° — D° mixing differs by removing the discrete 
ambiguity in 12 — > <f>\2 + 7T. 



II. FORMALISM 



can, in general, be written as 



Af = 
A l = 
= 

Af = 



4e+ i */[l+r / e^+*')], 
^[l + r/e^-'W], 



A T f e 



A T e i(A f -tf )[1 + r _ ei(S7 -^ (4) 



where and Aj are the magnitudes of the dominant 

SM tree-level contributions. The ratios 77 and r— are the 
relative magnitudes of subleading contributions contain- 
ing new weak phases (they could arise from new physics, 
or from SM amplitudes with suppressed CKM structure). 
<^J, (fiL : (f>f , and cf>j are weak (CP violating) phases which 

appear with opposite signs in CP conjugate amplitudes, 
and A/, 5f, and Sj are strong (CP conserving) phases 
which appear with the same signs in CP conjugate am- 
plitudes. All of the quantities entering Eq. ((4|, except 
the weak phases, are understood to be phase space de- 
pendent for 3-body and higher final states. 

In the case of decays to CP eigenstates, Af = 0(tt) for 
CP even (odd) final states. Eq. ((4|) therefore reduces to 



Af = 
Af = 



A T r ^ 



'If 



f [I + r/e' 

f[l + 
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where rfj p = + (—) for CP even (odd) final states. The 
time-dependent CP asymmetries depend on the universal 
quantity 



x _ qA f 



Vf 



i<t> 



(G) 



We begin with a summary of the formalism for neu- 
tral meson mixing and decays [l|, Q • The neutral meson 
mass eigenstates are linear combinations of the strong 
interaction eigenstates \M°) and \M°), 



|Mi, 2 ) = P \M°)±q\M°) 



(1) 



where \q\ 2 + \p\ 2 = 1. We define the mass and width 
differences as 



m 2 ~ m i 



y 



2T 



(2) 



where r/ of Eq. is neglected in the equality, and 
<f) is the relative weak phase between the mixing and 
decay amplitudes. Examples of decays to CP eigen- 
states include D° -> K + K~. tt+tt", K s tt , and B s -> 

J/*<j>, Di* )+ D^-. 

In the "pure-penguin" decay B s —* (fxf>, Aj is the mag- 
nitude of the Standard Model penguin amplitude. Ne- 
glecting r/, the weak phase <j) in Eq. ((SJ is the same as in 
the tree-level B s examples above, up to a small correction 
Scj> = 2Im(V: s V ttb /V*V cb ). 

For final states which are not CP eigenstates, the time- 
dependent CP asymmetries depend on 



where T = (Ti + T 2 )/2 is the average width. 

The decay amplitudes of the neutral mesons M° and 

M° to CP conjugate final state / and / are denoted as 



A f = (f\H\M°) , Af = (/|«|M°) , 
7 = (7\n\M°), Xj={f\H\M°), 



A 



(3) 



where TL is the weak interaction effective Hamiltonian. 
The decay amplitudes for the tree-level dominated decays 



A/ 



At = 



qAf_ 
pA f 
qAj 
pAj 



R1 



J 



-A f ) 



(7) 



where 77 and rj of Eq. (0J are neglected in the equalities, 
and thus Rj 1 = R f = Aj/Aj. Examples are D° -> 
K*% ± ,K*K and B s -> DfD* T . The weak phase (j) is 
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the same in Eqs. ((6|) and (O for the D° decays (up to 
negligible corrections < \(V uh V cb )/(V us V cs )\ ~ 10~ 3 ) and 
tree-le vel B s decays listed above. 

The M° — M° transition amplitudes are 



(M°\H\M°) = M 12 --r 12 . 
(W\H\M°) = M* 2 - l -T* 12 . 



(8) 



where H is the 2x2 effective Hamiltonian governing 
neutral meson mixing. We define the mixing parameters 

x 12 = 2\M 12 \/T, y 12 = |r 12 |/r, <j> 12 = arg(M 12 /r 12 ). (9) 

The notation x\2, yi2 is borrowed from [gj. 0i 2 is a CP vi- 
olating weak phase which is responsible for CP violation 
in mixing 7^ 1). Solving the eigenvalue problem 

yields 

(x-iy) 2 = x\ 2 - y\ 2 - i2x i2 y 12 cos <p 12 , or 

x 12 - Vl2 . X V = X \2V\2 COS 012 , (10) 



2 2 

x -y 



and 



9 _ -r(z - iy) 
P ~ 2(Mi 2 



-2(jgg - |T|g) 
T(x - iy) 



(11) 



The phase transformation \M°) -> e lS |M°), |M°) -> 
e _ie |M°) has no physical effects, due to conservation of 
Strangeness, Charm, or Beauty number by the strong in- 
teractions. Indeed, it is easily seen that </>i 2 , A/, Ay, x, 
and y (which are related to, or are themselves observ- 
ables) are invariant under these phase redefinitions [l|. 
Furthermore, the mass eigenstates are rotated by a com- 
mon phase factor. 

One is free to identify M 2 or Mi with either the short- 
lived meson (Ms) or the heavier meson (Mh), by re- 
defining q — > —q. This is equivalent to choosing a sign- 
convention for y, which in turn fixes the sign of x via 
sign(cos 0i 2 ), or vice- versa. Note that changing the sign- 
convention for y (or x) takes A/ — > —A/, or equivalently, 
<f> — > (/> + 7r. However, the combinations y A/ and x A/, or 
?/ cos 0, y sin and x cos 0, a; sin <j> are sign-convention in- 
dependent, which is seen explicitly from Eq. (fTTj) . Thus, 
they are candidates for being related to physical observ- 
ablcs. 

Examples of CP conserving observables are 
sign(y cos </>) and sign(x cos (f>) . In the limit of small 
or no CP violation, respectively: (i) Ms would be 
approximately or exactly CP-even if and only if 
sign(ycos0) = +1, and (ii) Mh would be approximately 
or exactly CP-even if and only if sign(xcos0) = +1. 
This is seen from Eqs. ([T]), (|3|), and (|6|) by requiring 
that CP-even (M + ) and CP-odd (M_) states decay into 
CP-even and CP-odd final states, respectively. In fact, 
in the D° system in the limit of CP conservation, the 
observable yep, defined in Eq. (f3"2"| . is equivalent to [|[ 



2/cp 



r(£> + )-r(£>_) 

r(r>+) + r(£>_) 



The world average is @, 

y C p = (1.07 ±0.26)%. 



(13) 



Taking into account that \q/p\ ~ 1 and |sin0| <c 1, see 
Eq. (j3"3"|) , one finds that yep ~ y cos <fi to very good 
approximation [4J, thus explicitly realizing (i) above. 

An alternative choice employed by the PDG [13] and 
HFAG L 4] collaborations for the A" and D° systems, is 
to identify M 2 with the would-be CP-even state in the 
limit of no CP violation. This amounts to choosing a 
convention for cj>, i.e., <j> ~ rather than (f> rj 7T. Given 
that in both systems the approximately CP-even state 
is Ms, this choice is equivalent to the sign-convention 
y>0. 

If M 2 were identified with Ms (y > 0), Eq. (10]) would 
give 

x = sign(cos0i 2 ) x 

x\ 2 - y\ 2 + yj (x\ 2 + y 2 12 ) 2 - ±x\ 2 y\ 2 sin 2 j> V2 j . ( 1.4 ) 

y 



Vxi ~ x ' 2 i2 + \l (aaa + vh) 2 ~ ^liVli ^ <fat\ • (15) 



If, instead, M 2 were identified with Mji (x > 0) , then the 
factor sign(cos0i 2 ) would be moved to the equation for 
y, with appropriate modifications for the choices y < 
or i < 0. These equations relate the the neutral me- 
son mass and width differences to the underlying mixing 
parameters xi 2 , y%2, and 0i 2 . 



III. THE CP ASYMMETRIES 

CP violation in pure mixing corresponds to \q/p\ =/= 1. 
It can be measured via the "wrong-sign" semileptonic CP 
asymmetry, 



(12) 



r(M°(t) -► t~x) - r(M°(t) -> e+x) 

T(M°(t) -> t-X) + T(M°(t) -> e+X) ' 
= (|^| 4 -l)/(|g/p| 4 + l). (16) 

In the limit \\q/p\ — 1 <C 1, which holds to good approx- 
imation for all four meson systems, ash — 2(\q/p\ — 1). 

The D° time-dependent decay rates into a final state 
/ can be written as (see, for example, Q) 

r(D°(t)-/) = \e- T \A f \ 2 {{l + \\ f \ 2 ) cosher) 
+ (1 - |A/| 2 )cos(a:T) +2Re(A/) 
X sinh(yr) - 2Im(A/)sin(a;T)|, (17) 

r(Dff(t)-/) - ie- T |3 / | 2 {(l + |A7 1 | 2 )coBh(yr) 
+ (1 - | AJ 1 1 2 ) cos(a;T) + 2Rc(Aj 1 ) 
x sinh(yT) - 2Im(Aj 1 ) sin(xr)| , (18) 
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where r = Tjjt. 

For D° decays to CP eigenstates the above expressions 
yield, to good approximation, purely exponential forms 
due to the small values of x and y, 



r(D°(i) ->/) cx expI-IV^ t], 
T(W(t)^f) cx exp[-f^ / t]. 

The decay rate parameters are Q 



(19) 



r D ^/ = T D [l + r]J \q / p\ (y cos (j)-x sin (/>)], 
f w ^ f = T D [1 + V f p \p/q\(y cos (I, + x sin (/>)}, (20) 

where <f> is defined in Eq. ©, and rj has been neglected. 
Note that Eq. (|20|) applies to singly Cabibbo suppressed 
(SCS) 2-body decays (e.g., D° -> K+K~ , tt+tt"), and 
to 2-body decays in which both Cabibbo favored (CF) 
and doubly Cabibbo suppressed (DCS) amplitudes con- 
tribute (e.g., D° — ► K s ir a ). (In the case of decays to CP 
eigenstates which are resonances or multi-body states, 
Eq. (f2TJ)) is valid when ignoring the interference of these 
amplitudes with other amplitudes in phase space, see be- 
low.) One defines the CP violating combination (or life- 
time CP asymmetry), 



AY f = D = a m + a\ (21) 



where 



2r D 



-^ p |cos< 



Vf 2 Sm< 
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q 




p 


) 




p 




q 





( 


q 


+ 


p 




p 




q 



(22) 



a m and a % are the contributions due to CPVMIX (\q/p\ ^ 
1) and CP VINT (sin<j> ^ 0), respectively, and are uni- 
versal quantities. Note that they are independent of sign 
convention for x or y. Subleading, non-universal correc- 
tions to T D a^j, Tjyt^j; due to rj ^ are discussed in 
Section VI. 

In SCS D° decays to non-CP eigenstates (e.g., D° — > 
K*K), the final states are essentially resonances or multi- 
body states. The time-dependence of the decays is again 
exponential, to good approximation, and is independent 
of phase space if the interference of these amplitudes with 
other amplitudes is ignored. In general, in decays to 
resonances, or multi-body decays, the exponential decay 
rate parameters depend on phase space (e.g., for 3-body 
decays, the location in the Dalitz plot) and give two CP 
violating combinations (llj |. 



AYf EE 



2r; 



= a T + a), 



f 2Tn f f 



(23) 



where, neglecting rf and ry, 

p y 'f 
—Rf — r cos < 



£ 
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P 
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«/ = R fY sli " 
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+ 


p 




p 







(24) 



(for aj, replace /->/), 



x'j = x cos Af + y sin Af, y'f = y cos Af — x sin A f , 
x j = x cos Af — y sin A/ , y'j—y cos A/ + x sin A f , 

(25) 

and 0, Af are defined in Eq. (0. In SCS decays one 
expects Rf = 0(1), implying that the CP asymmetries 
for non-CP eigenstates should be of same order as for CP 
eigenstates. 

The quantities af, a l f , , aL are not universal for 

non-CP eigenstate final states, due to the presence of 
strong phases. However, the latter can be determined, 
e.g., for 3-body decays, from Dalitz plot analyses. For 
example, in the simple case of a single resonance, K*K, 
in the Dalitz plot, Ak-k can be determined from the 
interference region of K + *K~ with K*~ K + [Hj]. Conse- 
quently, x, y, \q/p\, and <f> can be determined (up to dis- 
crete ambiguities) in Dalitz plot analyses of final states 
such as D° -► K 8 i^Ti-T and D° -> tt-tt+tt [ll|. 

In the case of CF and DCS decays to non-CP eigen- 
states, the time-dependence for D° decays to the "wrong- 
sign" (WS) final states D°(t) -> / and D°(t) -> / is ex- 
panded to quadratic order in r inside the curly brackets 
of Eqs. (fTT|) and (JTSJl, due to the small values of tan 2 9 C , 
x, and y (Aj is chosen to be the DCS amplitude, e.g., 

D° — » K + ir~ or / = K + ir~). The result can, in general, 
be written as (we adopt a notation similar to the one 
used in the experimental analysis of D° — > K^tt^ [IF 



Y\D\t) - /] = e^|A / | 2 



(R+f + R+y'+r 



T[D°(t) - /] = e^\Aj\ 2 x 



(RJ) 2 + RJ y'- t 



(26) 



Qr'+) 2 + (y'+f r2 



(27) 



(x'-f + (y'- f r2 



where R't = \Aj/Af\ and Rf = \Af/Aj\ are the magni- 
tudes of the DCS to CF amplitude ratios for D and 
D° decays. Neglecting r/ and rj (until Section VI), 

Rf = RJ = Rf = \Af/Af\ = Aj I Aj , as in Eq. 0, 
so that 



V* = 



x'± = 



q 




p 


p 




(i 


q 




p 


p 




q 



x (x'j sin 4>Tyj cos 0), 

x (a4 cos0 ± y'j sin0), (28) 
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where x'j and y'j have been defined in Eq. (I25|) . In addi- 
tion, 



(29) 



allowing to be expressed solely in terms of {x'^) 2 + 
iy'^) 2 above. The expressions given in Refs. [H and [l3| 
for y /=t and x ;± differ from those in Eq. (j28|) due to choice 
of convention, and are recovered by substituting Af — > 
— A_f + 7r in and y^. The time dependence for D° 

decays to the "right-sign" (RS) final states is, to good 
approximation, exponential and given by 



T[D°{t) 



f] = 
7] = 



2 



~\ A 7 



(30) 



Thus, the decay rate parameter is r d°^k- tt+ = Td- 

A fit to the time-dependence in Eqs. (|2"B")) and (j2"T|) 
yields measurements of , y'^, and x ,:iz , which can be 
used to determine or constrain 1 — \q/p\ and </>, as carried 
out in [13] for D° -> K ± n^. Note that the CP violating 
quantity (y 



satisfies 



Rj(y' + -y'-) = AY- ] 



f 



(31) 



where Rj = R^ is the magnitude of the CF to DCS 
amplitude ratio (for rf = rj = 0), see Eqs. (O and 
(HUMUS)- Finally, the contributions of CPVMIX and 
CPVINT in D° decays to RS final states are relatively 
suppressed by tan 4 8 C , and are therefore not considered. 

An important CP conserving quantity yep , mentioned 
in Section II, can be defined in terms of the decay rate 
parameters T D o^f CP (for SCS decays to CP eigenstates) 

and r.D ^.R"-7r+ j 

ycp-Vf -» 1- (32) 

The expressions for the decay rate parameters given 
above (in the rj = limit) imply [i[ 



Vcp = ^ cos < 



£ 




P 




£ 




P 
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)-§sm( 
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P 




Q 




P 




Q 





.(33) 



The time-integrated CP asymmetry for D° decays to 
CP eigenstates (SCS and CF/DCS) is defined as 

af ^(D^.f)-T(^f)_ (34) 
T(D° -> /) +T(D° -> /) 

Expanding to leading order in x, y, rj yields [ll[ 

a f = aj + a m + a 1 , (35) 

where a m and a 1 are given in Eq. ([22]) . and 

a d f — 2r f sin (j) f sin S f (36) 



is the (non-universal) direct CP violation contribution. 
The time-dependent CP asymmetry (AYf) and the time- 
integrated CP asymmetry (a/) are equal if there is no 
direct CP violation. 

For SCS D° decays to non-CP eigenstates there are 
two time-integrated CP asymmetries to consider, 



r( d° -> /) - r(W -> /) 
_ t(d° - /) - r(^ -» /) 



°7 = 



(37) 



r(£° -f /) + r(£° - /) 

Expanding to leading order in x, y, rf, ry yields [ll[ 

a/ = aj + af + a), aj = aj + a'j + Oj, (38) 
where a™, Oj and Of, aL are given in Eq. and 

a j = 2r/ sin / sin 5/ , a 4 = 2rj sin cj>j sin <5j (39) 

are the direct CP violation contributions. Again, if there 
are no new weak phases in decay, the time-dependent and 
time-integrated CP asymmetries are equal, i.e., AYf = 
a f and AYj = aj. 

In the case of CF/DCS decays to non-CP eigenstates, 
and in our convention for RS and WS final states, the 
definitions of a/ and aj in Eq. (f3"T|) correspond to the RS 
and WS time-integrated CP asymmetries, respectively 
(e.g., a K -K+ an d a K+iT- f° r D° — > K^^). To leading 
order in x, y, rf, and rj they are given by 



l% aj = Rj(y ,+ -y 



(40) 



where the RS (ai) and WS ((£) direct CP asymmetries 
are as in Eq. (|39j) . 

The timc-dcpcndcnt CP asymmetry for B s decay to a 
CP eigenstate, to leading order in rj and for \q/p\ = 1 
(the HFAG average is \q/p\ = 1.002 ± 0.005 (g), takes 
the simple form [l| 

r(B,(t)^/)-r(g(t)^/) 
r(B a (t)->/) + r(s, (*)->/) 

= 5/sin(|a;|rt) -C/cos(|x|rt), (41) 

where 

5/ = 77? p sign(ir) sin0 , C/ = 2r/ sin0/ sin 5/ (42) 

are the contributions due to interference between mixing 
and decay, and direct CP violation, respectively. The 
factor sign(cc) in Sf originates from the time-dependence 
of the decay rates, via sin(xrt) = sign(x) sin(|x|r£), and 
insures that Sf is independent of sign convention. 

For B s decays to non-CP eigenstates there are two 
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time-dependent CP asymmetries to consider, 
T(B s (t) -> f) - TjBJf) -> 7) 

r(B a (f)-»/) + r(BJ(t)->/) 

= S f sin(|x|rt) - C f cos(|x|rt) , 

r(B s (t) -» J) - r(TO -> /) 
r(s a (t)-/) + r(B.(t)-/) 

= 5jsin(|x|rt) -Cjcos(|x|rt), (43) 

where (again to leading order in r*y, and for \q/p\ = 1), 

S/ = Sj — 2 sign(x) sin(0) cos(5/) r//(l + r^) , 

Cj = 2rj sin 0j sin Sj . (44) 



Cf = 2rr sin (if sin(5f 
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The equality between 5/ and Sj holds, up to negligible 
corrections of 0(\q/p\ — 1). 



IV. RELATING THE INDIRECT CP 
ASYMMETRIES 

In general, we are i nter ested in decays to final states 
common to M° and A/ , whose leading contributions 
to Ti2 are proportional to the dominant CKM struc- 
ture entering this quantity, i.e., (V CS V* S ) 2 for the D° and 
(V cb V* s ) 2 for the B s . All of the examples we have men- 
tioned previously are in this class. In this section we as- 
sume that there are no subleading amplitudes with new 
weak phases in these decays [rf = rj = in Eq. Q], and 
we neglect CKM suppressed contributions to T^. The 
following relations are then satisfied: 



is related to CP violation in mixing. Note that Eq. I|48[). 
which also appears in Q , relates CP VMIX to the under- 
lying mixing parameters x%%, 2/12, and 0i2. 
Multiplying (see Eq. ([TT]) ) 



M- 



12 



2 1 12 



Ml2 



(51) 



on the l.h.s. by (Af/Af) 2 for decays to CP eigenstates to 
obtain \ 2 (or by (AfAj)/(AfAj) for decays to non-CP 
eigenstates to obtain A/ Ay), and on the r.h.s. by T^/r^ 
yields 



tan 20 



sin 2(f> - 



sin 2(f>i 



COS 2012 + Vlz/xii 



2A m xy 

y 2 + A 2 m x 



cos 2(f) 



y 



A 2 x 2 



2 + A 2 x 2 



y 



(52) 
(53) 



The first relation is incorporated into the fit of X12, 
2/i2, and 0i2 using the D° — D° mixing data. The 
last two relations are obtained by eliminating the de- 
pendence of sin 2(f> and cos 20 on X12, 2/12, and 0i2, us- 
ing Eqs. (fTT)|) . (|47H49| . Finally, a trigonometric identity 
yields 

tan0 = -A m x/y . (54) 

This expression also appears in Q ■ It relates CP VMIX to 
CPVINT, model-independently, in decay modes in which 
there are no new weak phases, and is independent of sign 
convention for x or y. In the limit \\q/p\ — 1| << 1, which 
holds to very good approximation for all four meson sys- 
tem, we obtain 



A f A } 
A)A S 



6l 

A f 



and 



A f A *f 



A f A J 



r* 2 A'fAf + AfA 



f 



A f Aj 



(45) 



(46) 



for CP-eigenstate and non-CP eigenstate final states, re- 
spectively. CKM suppressed contributions to Ti2 and to 
Tf,rj within the SM yield corrections to these relations 
of 0(\(V cb V ub )/(V cs V us )\) « 6 • 1(T 4 for D° decays, and 
of 0{\(V ub V us )/{V cb V cs )\) « 0.02 for B s decays [see Eq. 

GEM. 

The following formulae, obtained from Eqs. 
and (jlip. will be useful: 



\q/p\ 2 {x 2 +y< 

4 



x u + 2/i2 + 2x122/12 sin 0i2 



Vl2 



x 12 



2/12 



2xi22/i2 sin 012 



x i2 + 2/i2 _ 2xi22/i2 sin 012 



V 2 + A 2 m x 2 
1 A 2 



A 2 m y 2 



b 12 



1 — A 2 



(47) 
(48) 

(49) 



where 



tan 0=1 — 



x _ a SL x 
V 2j/ 



(55) 



As discussed in [5j, this relation gives an excellent de- 
scription of the data in the neutral kaon system. 

It is straightforward to relate AYf and the semilep- 
tonic CP asymmetry using Eq. (|54[) . after expanding to 
first order in \q/p\ — 1. In the case of D° decays, the same 
relations also apply to the time-integrated CP asymme- 



tries (for rf 
one obtains 



0). For decays to CP eigenstates, 



a cp fl SL y 2 +x 2 

AY f = a f = -ycos^Tjy — 2 — . (56) 

^ 2/ 

We know from experiment that the level of CP viola- 
tion in the D° system is small and that the short-lived 
meson is approximately CP-even, implying |cos0| w 1, 
sign(j/cos0) = +1 (as in the Standard Model) and, to 
good approximation, 



Av C p a shy 2 + x 2 

*Tf = «/ = ->/ - T 



(57) 



which is independent of sign convention for x or y. Sim- 
ilarly, we obtain 



An 



{\q/p\ 2 -l)/{\q/p\- 



(50) 



2/cp = y/ cos0 = \y\ , 



(58) 
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up to corrections of order sin or a 



For SCS D° decays to non-CP eigenstates, one obtains 



R f AY 7 = AYf/Rf = R f aj = a f /R f 



cos A 



asL y 



\y\ 



(59) 



Confirmation of the relation between the hadronic CP 
asymmetries in the first line of Eq. (|59|) does not require 
knowledge of A, . In terms of the CP-averaged branching 
ratios for D° — > / and D° — > / decays, it is simply given 
by 



ay-. 



/ 



AY, 



f 



«./■ 



Br(£>° -> /) 
Br(£>° -> /) 



This relation follows non-trivially from Eq. ([5 

,C'P\ A V. / D . I ^ „CP I T> .\ could 



(60) 



RtAYf 



f f 

(or Rf and AYf/Rf (or a^ p /Rf) could, in princi- 

ple, differ by 0(1) given that y ~ x, that sin A/ could be 
large, and that (|g/p| — 1) ~ sin0 is allowed. 

In the case of CF/DCS D° decays to non-CP eigen- 
states, Eqs. (HU), (gO}, and (J59J) imply that the time- 
dependent, time-integrated, and scmilcptonic CP asym- 
metries are related as (recall / = K~ir + for D — ► Kir in 
our convention) 



= Rftij = — cos A / 



asLp ) 



(61) 



The strong phase Ak-k for D° — » K ± tt t decays can be 
precisely measured by the BES-III Collaboration at the 
^(3770) charm threshold. 

For B s decays to CP eigenstates, the time-dependent 
and semileptonic CP asymmetries are related as 



,1/2 



-rjf p sign(y cos < 



)a SL k/y|, (62) 



which is independent of sign convention for x or y (\x\ 
follows from sign(x) in Sf). At this point, we elaborate 
on the determination of sign(y cos 0) in [5j . Starting with 
Eq. (jlip . the last relation in Eq. ((TDJ) , and taking y <C x 
(using the HFAG averages [9j, the central value for y/x 
is w 0.003), we obtain 

sign(ycos0) = sign(|Mi 2 | 2 Re[rja4/M/] + 

Im[Mi 2 rt 2 ]Im[M 1 * 2 A / /A / ]). (63) 

The ratio of second to first terms above is given by 
sin sin 0/ cos (0i2 +0). However, for y <C x, 0i 2 = 
0mod(7r), see Eq. (|52j) . implying that the magnitude of 
the ratio is less than 1. Thus, Eq. (f6"3")) simplifies, and 



sign(?/cos0) = sign(Re[T* 12 A f /A f }) . 



(64) 



Given that the impact of new physics on the r.h.s. would 
in general be subleading, we conclude that sign(y cos 0) = 
sign(y cos 0)sm = +1> and that 



25//(l 



s 2 s y' 2 



Vf \x/y\a S h 



(65) 



x[%] 


v l%] Wp\ 


[rad] 


1.00 ±0.25 


0.77 ±0.18 0.94 ±0.14 


-0.046 ± 0.093 


1.00 ±0.25 


0.76 ±0.18 0.86 ±0.16 


-0.15 ±0.13 


6 K n [rad] 


Skwtt [rad] R D [%] 


A D [%] 


0.40 ±0.19 


0.20 ±0.37 0.336 ±0.008 





0.39 ±0.18 


0.20 ±0.37 0.336 ±0.009 


-2.1 ±2.4 



TABLE I: HFAG outputs for Ad = (first row) and Ad ^ 
(second row) 



in the absence of new weak phases in decay, as in [f| . 
For decays to non-CP eigenstates, 



Sj = -k \x/y\ a SL /2 , 



(66) 



where 



(AR) cos z A f - Sj {R) ± l) 2 ) 1/2 /(Rj ± 1) , (67) 



and, as usual, Rf = \A f /A f \ = Al./A T f . The (near) 

equality of the two time-dependent CP asymmetries, al- 
ready noted in Eq. (|44|) . is a trivial consequence of y <C x 
and \ \q/p\ — 1| <C 1, unlike in D° decays. 



V. FIT RESULTS FOR D° - D MIXING 

The current D° — D mixing and CP violation fit re- 
sults reported by the Heavy Flavor Averaging Group 
(HFAG) Q can be expressed in terms of the four uni- 
versal parameters (x, y, \q/p\, and 0), two strong phases 
(8 Kir and 8Kirir)i the CP averaged ratio of wrong-sign to 
right-sign D° — > if ± 7r T decay rates (Rd), and the corre- 
sponding direct CP violation parameter (Ad)- In terms 
of our notation for CF/DCS decays, 



Rd 



An = 



(Rj) 2 



(R+) 2 ~ (RJ) 2 



(RI 



(RJ) 2 



(68) 



with / = K~it + . The four universal parameters are ex- 
tracted from fits to the time-dependent decay rates for 
D° — » K + K~, it + it~, K^tt^ , Kirn, and the semileptonic 
decay rates The HFAG fit only allows for new weak 
phases in the D° — > K ±r K T amplitudes, via Ad ^ 0. New 

weak phases in decay and their impact on the D° — D 
mixing and CP violation fit are discussed in more detail 
in Sec. VI. 

In general, x, y, and \q/p\ can be expressed in 
terms of the mixing parameters Xi 2 , j/i 2 , 0i 2 , see 
Eqs. HU), (fl5|) . (|4"8"|) . In the absence of new weak phases 
in decay, the same is true for 0, see Eq. (|52|) . Using 
these four equations (recall that Eqs. (TT4")) and (fTS)) cor- 
respond to the HFAG convention, which identifies M 2 
with the approximately CP-even state) x, y, \q/p\, and 
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Parameter 


A D = 


Eq. ([52l removed 


A D 


X 12 [%} 


1.00 ±0.25 


1.00 ±0.25 


1.02 ±0.24 


yi2 {%] 


0.77 ±0.18 


0.78 ±0.18 


0.75 ±0.18 


0i2 [rad] 


0.02 ±0.08 


-0.12 ±0.22 


0.07 ±0.08 



TABLE II: Results for the mixing parameters at la, see Sec. V 
(VI) for A D = (A D 0). 



are determined by the mixing parameters X12, j/12, 012- 
Ranges for these underlying parameters can be extracted 
directly from experimental data under the assumption 
that An = 0; where HFAG currently reports seven pa- 
rameters, only six would be reported. 

For this work, we adopt a simpler strategy for extract- 
ing values of xu, j/12, and </>i2: we take the HFAG fit 
results (for the Ap = case) for x, y, \q/p\, <fi, ok-k, 
ok-k-k, R-D, shown in Table HI and minimize 



A" 



(69) 



where is the difference between the HFAG value for 
the i th parameter and the fitted value predicted using 
the equations which relate x, y, \q/p\, and </> to £12, 2/12, 
and 4>i2] the weight matrix Wij is the inverse of the full 
error matrix for the values reported by HFAG, including 
the correlation coefficients [lj]. The fitted values for &k-k, 
Sktttt, and Rd are very close to the HFAG values; they 
change only due to (small) off-diagonal elements in Wij . 
The HFAG parameters used as input are taken from a fit 
with x 2 = 24.9 for 21 degrees of freedom (28 experimental 
results minus 7 parameters). The value of \ 2 in our fit 
is 0.2. In effect, the overall x 2 increases slightly as one 
degree of freedom is restored to the mix of measurements 
and the parameters to be extracted. 

The fitted values of X12, yi2, and 4>\2 are listed in the 
second column of Table HT1 In particular, we obtain, 



[rad] = 0.02 ±0.08 (la). 



(70) 



Our results for 212 and j/12 are very close to the fitted val- 
ues for x and y in Table HI as would be expected for small 
</>i2, see Eqs. (|14p . (fT5|) . A bound equivalent to a precision 
on (f>f 2 of ±0.18 (la), which assumes no correlations be- 
tween the ex per imental measurements, has recently been 
obtained in hj. The HFAG error matrix corresponds 
to parabolic errors, and thus our two sigma and higher 
CL intervals are simple multiples of our la CL interval. 
However, a preliminary HFAG fit [T(J| to the data, which 
uses Eqs. (fl"I)l , (fT5|) . P5|) , and (j5l2|) , as discussed above, 
indicates that the errors on 0i2 are non-parabolic (and 
thus we do not list higher-CL intervals). Therefore, our 
fit result for ^12 is only approximate. The preliminary 
HFAG la and 95% CL intervals for non-parabolic errors 
are 



. [rad] = 0.02 



+0.06 
0.13 



(la) 



The former is similar to our result using parabolic er- 
rors. The HFAG fit results for parabolic errors are in 
agreement with ours. 

The impact of the relation between (j) and <f>\2 on the 
precision with which 4>\2 is constrained is seen by re- 
peating the fit for the Ao = case, but with Eq. (j52"|) 
removed. In this case 4> is treated as an independent 
parameter which is trivially fit. The result is reported 
in the third column of Table [Til We observe that the 
error on ^12 increases by roughly a factor of three, and 
thus conclude that the relationship between CPVMIX 
and CPVINT provides a powerful constraint on the al- 
lowed mag nitude of CP violation in D° - D° mixing. 

Finally, to understand the implications of the bound 
on 4>i 2 for model building, we separate M12 into its SM 
and new physics parts, 



M12 = Mf-fe^™ 



(72) 



SM 



where only the difference of the weak phases <f)^f — y> M 
is physical. We continue to assume that there are no new 
weak phases in decay, and identify Ti2 with its SM value. 
The definition of 0i2 then yields 



sin 0f 2 



A/ NP 



Mis 



(73) 



where IM12I follows from the fitted value of xi2- In the 
usual phase convention in which Aff 2 M is real (^f/ 1 = 0), 
the above bounds on <^ thus imply that 

Im(M^) 



|Mi 2 

for parabolic errors, and 
Im(M^ F ) 



e [-0.06, ±0.10] (la) 



(74) 



|M 



12 



e [-0.11, ±0.08] (la), 
E [-0.30, ±0.30] (95% CL). (75) 



G [-0.30, ±0.30] (95% CL). (71) 



for the (preliminary) non-parabolic HFAG errors. As 
shown in the next section, these bounds can not be sub- 
stantially altered if we allow for new weak phases in de- 
cay. 



VI. NEW WEAK PHASES IN DECAY 
A. General considerations 

In this section we discuss how the relations between 
CPVMIX and CPVINT are modified by new weak phases 
from subleading decay amplitudes (originating from new 
physics, or CKM suppressed SM amplitudes). We be- 
gin with the resulting shifts in arg(A/), arg(Ay), and 
arg(r:f 2 /ri2). Expressions relating arg(A/) and arg(Ay) 
to 1 — \q/p\, as well as to (f>i2, which depend on these 
shifts, are obtained, replacing the previous expressions 
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involving <f>. In turn, new relations between the time- 
dependent and semileptonic CP asymmetries are derived 
for D° and B s decays. Direct CP violation bounds are 
used to constrain deviations from the r/ = r j = case, 
and the la intervals for X12, 1/12, and 0i2 from an ap- 
propriately modified fit to the D° — D° mixing data are 
presented. 

The argument (j>\, = arg(— A/) for a decay to a CP 
eigenstate in Eq. ([6]) is shifted, to first order in r/, as 

a f = <t> + $<j>\ f , S<fi\f = -2r/ cos (5/ sin 0/ . (76) 

For non-CP eigenstates, the arguments <f>\, = arg(— A/) 
and (f>\f = arg(— Aj) in Eq. (7J) are shifted by 

8<f>x f = -rfsin(6f + <f>f) + rjsm(§j-<f>j), 
8<j>xj = -rjsin(Sj + 0j) + 77 sin(<5/ - <j) f ). (77) 

The new contribution to Arg(ri2/r| 2 ) is defined as 



arg 



12 



r 
1 ii 



1 12 
1 12 



= 21m 



££12 

1 12 



(78) 



to leading order in <5Pi2 = Ti2 — r° 2 , where r° 2 is the 
leading SM contribution to Ti2 proportional to (V CS V* S ) 2 
for the D° and (V cb V* s ) 2 for the B s . Note that <50 r is 
phase redefinition invariant and is an observable, unlike 
arg(ri 2 ). 

84>v receives contributions from CKM suppressed cor- 
rections to Ti2 within the SM, and from subleading de- 
cay amplitudes (rf,rj ^ 0). (We note that a recent 
analysis of Ti 2 in the D° system [TtJ indicates that the 
CKM suppressed corrections to Sfir could be enhanced 
from 0(\V cb V ub /V cs V us \) in the SM to 0(0.01) in mod- 
els with a fourth family.) The contribution to S(f>r from 
subleading decay amplitudes (£0f)j expressed as a sum 
over exclusive final states, and to leading order in r/, ry, 
is given by 



£ P ( A / f r S co s Sf sin 0/ 



(79) 



A^Aj(rf cos(A/ — 5/) sin0/ + rycos(A/ + (5y) sin0y), 

where the sums are over CP and non-CP eigenstates. 
We learn that 8(f>^ is of 0(4f/ sin 0/), roughly weighted 
by the fraction of that is attributed to the affected 
decay amplitudes within the SM. f/ is the "generic" size 
of Tf and rj in these amplitudes. The same qualitative 
conclusion can also be reached via the OPE treatment 
for Ti2, in the case of the heavier Bd and B s mesons. 
The relation between and 0i2 in Eq. (|52)l is replaced 

by 



tan(20 A/ - 25(j)\ f 







sin 20 



12 



cos 20i 2 + y 2 2 /xf 



(80) 



for decay to a CP eigenstate. The argument on the 
l.h.s. is simply 20 + <50r, which takes into account the 
shift in arg(ri2/r^ 2 ) in Eq. (|45| . The relation between 
CPVMIX and CPVINT for decay to a CP eigenstate is 
now given, in terms of the observable <j)\ f , by 

tan(0 A/ - 54> X/ + Hv/2) = -A m x/y . (81) 

Expanding to lowest order in r/ and \q/p\ — 1 yields 



tan A/ 



1 



y 



cos^ 



(82) 



Corrections to the relations between the semileptonic 
and time-dependent CP asymmetries follow straightfor- 
wardly, see below. 

In the case of non-CP eigenstate final states, new re- 
lations which combine the effects of new weak phases in 
decays to CP conjugate pairs (thus removing the depen- 
dence on the strong phase A/) are obtained by substitut- 
ing 4>\ f — * (<fi\ f + 

<f>x T )/2 and «50 A/ (<50 A/ + Hx 7 )/2 
in Eqs. (f80l) - (f82"f . In practice, it may be more useful 
to consider their effect on each decay separately (see the 
discussion of D° — > K ± tt t decays below), yielding 



tan(20 A/ - 26<p Xf - 2A 



sin 2(f) 



12 



cos 20i2 + yi 2 / x 



12 



= tan(20A 7 -250A T + 2A/ + 50 r ) 



(83) 



for the dependence of the observables <j>\ f and <f)\- on 
012, and 



tan(0 A/ - 84>x f ~ A/ + 8<j> T /2) = -A rn x/y 

= ta,n{<t>x T -8<t>x T + &f+84> r /2) (84) 

CPVMIX and 



for D 



o 



V T "1 

for the modified relations between 
CPVINT. 

Approximate bounds on <50 A/ , 8<pXj, and <50r 
and B s decays can be obtained from direct CP violation 
measurements. It is instructive to compare them to the 
current experimental sensitivity to 8<f>x f , <50a^, and 8<fir 
in time-dependent (mixing-related) measurements. 



B. D° -D° mixing 

We need to consider new weak phases in singly 
Cabibbo suppressed (SCS) decays, and their combined 
effects in Cabibbo favored (CF) and doubly Cabibbo sup- 
pressed (DCS) decays. We begin with a discussion of the 
former. The HFAG average for AY/ obtained from 
the BaBar and Belle D° -> K+K~ and D° -> tt+tt" 
measurements 1181. is 



AY 



-0.123 ±0.248)%. 



The time integrated CP asymmetries for D° 



and D 



ii 



7r 7r are 

&K+K- 



9], 



(-0.16 ±0.23)%, 
(0.22 ±0.37)%. 



(85) 
K+K- 

(86) 
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The direct CP asymmetries are obtained by subtract- 
ing AYf from the time integrated CP asymmetries, see 
Eq. (J35J) , yielding 



4r+if- = (-0.04 ±0.34)%, 
<4 +7r _ = (0.34 ±0.45)%. 



(87) 



(Predictions for aj in the Standard Model suffer from 
large hadronic uncertainties spanning an order of mag- 
nitude or more, and could be as large as « 0.1%). Un- 
less the new physics has a very special structure, e.g., 
parity conserving [ijj], the results for aj (+K - and a^ +Tr _ 
give rough bounds on the direct CP asymmetries for all 
decays mediated by c — > u(ss, dd) transitions. Models 
which can easily produce direct CP asymmetries of this 
size or larger in SCS decays have been discussed in 

In general, strong phase differences enter as sin<5 in 
the direct CP asymmetries, and as cosi5 in 8<f>\ f , 8<f>\j, 
and 8<j>T- However, the relevant new physics (AC = 1) 
effective operators for SCS decays differ from the tree- 
level SM operators in their color and chirality structures 
(the QCD penguin operators, most notably the chromo- 
magnetic dipole operator, are relatively unconstrained by 
D° — D° mixing). Thus, strong phase suppression is not 



expected [11], implying that 
taking 



This justifies 



|^ A/ |,|<50 A _|,|^r| <1% 



for SCS decays, and similarly for the last term in Eq. ([82]) . 

The SCS decays enter the HFAG D° -W mixing fit 
via AYV (averaged over Tr + n~ and K + K~) and yep- In 
the case of decays to CP eigenstates a new weak phase 
would shift AYf, to lowest order in (1 — |g/p|) and rt, by 



8(AY f ) = -r,^(\y\a d f -\x\54>x f ) 



(89) 



This result follows by substituting (f> — > <j)\ f and \q/p\ — > 
\qAf /pAf \ in Eq. (|22[) . and expanding in small quan- 
tities. Note that the impact of r/sin0/ is suppressed 
by mixing (x,y ~ 10~ 2 ), unlike in aj which enters the 

time integrated CP asymmetry. With aj, 8<p\ f < 1% 
and x, y ~ 1%, we find 

|5(AY»| < 10" 4 , (90) 

which is less than a few percent of the experimental un- 
certainty, see Eq. (|85j) . The shift in yep due to new 
weak phases in decay must be even smaller relative to its 
experimental uncertainty, given in Eq. (|13[) . because its 
dependence on CP violating quantities must be quadratic 
(and still suppressed by x or y). 

The relation between ash and AY/ in Eq. ([ST]) for 
decays to CP eigenstates is modified, to lowest order in 
Tj and (\q/p\ — 1), as 



ash y 



2 

\x\{8<p\< 



\y\ a f 



Given that the approximately CP-even D° mass eigen- 
state is the shorter-lived and heavier one, we have sub- 
stituted ycoscj) — > |y| (as before) and xcos<j) — > \x\, and 
similarly below. Applied to SCS decays, the new physics 
correction is again < 10 -4 . The modified relations sat- 
isfed by ash, AY/, and AYj for SCS decays to non-CP 
eigenstates which replace Eq. (I5T)1) are 



AY f /Rf + RfAYj 
cos Af 



AYf/Rf - R f AYj 
sin Af 



-asL- 



V 



+ 



+ 



\v\ («/ 



\x\{6<j>xi +<50A T -<50r), (92) 
\y\{8(j>\ f + 8<p\ T - 8<p r ) 
\x\ (aj + af) , (93) 



1), and neglecting terms 
and are thus similarly 



to lowest order in r/, {\q/p\ 
of 0{rfX sin cj>), 0(r/ysin^ 
bounded. 

The most precisely measured CF and DCS time- 
integrated and direct CP asymmetries are near zero, 
with uncertainties of w 1% and w 5%, respectively. 
For example, the time integrated CP asymmetries for 
D° -> K-ir+n (CF) and D° -> K+tt-jt (DCS) are [| 

a f = (+0.16 ±0.89)%; / = /sTtt+tt , 
a f = (-1.4 ±5.2)%; f = K + n~n 



(94) 



The difference between the CF and DCS direct CP asym- 
metries in D° — * K ± tt t {Ap\ averaged over the BaBar 
and Belle measurements [lj, [20( , is 



= (0.4 ±3.5)% 



(95) 



and the global HFAG D° - D° mixing fit gives (-2.1 ± 
2.4)% (second row, Table I). Finally, the best CF D ± 
direct CP asymmetry bounds are for D + — > K~tt + tt + 
and D+ -> X s 7r+7r° ||, 



(-0.5 ± 1.0)%; 
(+0.3 ±0.9)%; 



/ 
/ 



(96) 



(tyr/2) 



(91) 



It is difficult to construct models with non-negligible new 
weak phases in CF and DCS decays [2l[ (one exam- 
ple is known [13]). Again, AC = 1 effective operators 
with different color and chirality structures than their 
SM counterparts would be important, and we can ex- 
pect significant strong phase differences. The direct and 
time-integrated CP violation bounds therefore imply 

\8<t> Xf \, \84> X7 \, \6<fr\ < 0(few percent) (97) 

(following our convention, take r/ and rj in 

Eqs. {32]) and (CB t0 correspond to CF and DCS 
new physics amplitudes, respectively). For complete- 
ness, we note that for CF and DCS decays to CP 
eigenstates (e.g., D° -> K s n°, p°K s ), \8<j> Xf \ < 0(1%), 
which is the approximate bound on CF direct CP 
violation (DCS contributions are suppressed by # 2 ). 
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New weak phases in CF and DCS transitions would 
enter the HFAG D° - 5° mixing fit via D° -> K^ir* 
and D° — » K s ir + ir~. For illustrative purposes, lets 
consider D Kir in more detail. The general form 
for the time-dependent amplitudes D°(i) — > i£T + 7r~ and 
Z>°(i) -> i^ _ 7r+ is the same as in Eqs. (HU) and (|2"7|). 
However, the corrected expressions for y and a; 1 * 1 , see 
Eq. (J2SJI, are given by (/ = #"77+) 



(* 






^/P 










(• 










A f p 




Ajq 



=£- x (sy sin^ =F 2/y cos^), 



x (a; j cos ± ± j/j sin <f> ) , 
(98) 

where (in terms of the direct CP asymmetry for the CF 
decays) , 



Aj/A f 



1 



d 
a f 



1 + 2r/ sin^/sin^/ , (99) 



and 



0a t + A/ 
*l - A 



(100) 



Corrections to the relation between the time-dependent 
CP asymmetry (y' + — y'~) = RfAYj and ash in Eq. |6T|) 
are easily obtained from Eqs. (f9"2")l and applied to 

CF/DCS decays. 



Measurements of y , a; and Rd, Ad [defined in Eq. 
([55)1 ] for D° — > Ktt have been reported by BaBar and 
Belle [HI, [2(| , also see @ . Averaging over the two exper- 
iments yields 



y' + - y'~ = (-0.19 ±0.64)% 



(101) 



for the time-dependent CP violation. The experimental 
error is an order of magnitude larger than the maximal 
allowed shift due to new weak phases in decay, of order x 
or y times the bound in Eq. (|9"7) . In addition, a fit for x'-j, 

yL and <ft has been carried out in the Belle analysis [l 

yielding 



(0.16 ±0.41) [rad]. 



(102) 



X \ of a few 



However, the fit uses the formulae for x'^ and j/ /=t in 
Eq. (|28p. thus neglecting the corrections in Eq. (f9"8")l . In 
particular, it assumes that <p + = <p~ — <f>. Fortunately, 
the reported error on is an order of magnitude larger 
than the upper bounds on \5<p\ f \ and 

percent, in 

negligible in Eq. (|9"9"]1. Thus, the use of Eq. Q28J) turns 
out to be a good approximation. 

The Belle Collaboration also fits for in a time- 
dependent Dalitz plot analysis for D° — > K s tt + it~ [23| . 
obtaining 



± . Moreover, a,j should be < 1%, hence 



(-0.24 ±0.32) [rad] . 



(103) 



Again, this analysis assumes that (f> + — $T — <f> (in 
general, <fi + and <\>~ would vary across the Dalitz plot). 
Again, the error on (f> is about an order of magnitude 
larger than the allowed shifts in ^ , implying that this 
is a good approximation. 

The outputs of the HFAG fit for An ^ (new weak 
phases in decay) listed in Table I have been obtained 
under the assumption that <p + — <\f~ = 4> in the time- 
dependent D° — » i4T ± 7r =F and D° — > K s n + n~ amplitudes. 
We have just seen that this is a good approximation. In 
addition, HFAG has not allowed for new weak phases in 
D° — * K + K~ and D° — * tt + it~. Again, this is a good 
approximation for SCS decays, given that the impact of 
new weak phases on AYj and yep would be negligible. 
Finally, modifications to the relation between 4> and 4>\2 
in Eq. ^ [see Eqs. (JHDJ) and ([55])] are smaller than the 
experimental sensitivity to <\> in CF/DCS decays by an 
order of magnitude, and in SCS decays by more than an 
order of magnitude. 

In view of the above considerations and in the case of 
new weak phases in decay, the mixing parameters X\i, 
yi2, and <\>\i can be obtained, to good approximation, 
along the lines of the fit carried out in Sec. V (for Ad = 
0). In particular, <p is once again expressed in terms of 
x\2, 2/12, and </>i 2 using Eq. (j52"| [as are x, y, and \q/p\, 
using Eqs. (fl"4"|) , (fT5)) . and l|48p]. However, now we take 
the HFAG fit results for Ad 0, see Table I, and add 
Ad to the sum over HFAG outputs in Eq. (f69|) for the 
X 2 function. The validity of this approximation reflects 
the suppression due to mixing (x or y) of the effects of 
new weak phases in decay on CPVINT [continued use of 
Eq. (|52|) ]. and the lack of such suppression in the direct 
CP asymmetries [use of the Ad ^ fit results]. The 
HFAG parameters used as input in this case are taken 
from a fit with \ 2 = 25.3 for 20 degrees of freedom (28 
experimental results minus 8 parameters). The value of 
X 2 in our fit is 1.3. Thus, as in the Ad = fit, the overall 
X 2 increases by a small amount as the number of degrees 
of freedom is increased by one. 

The fitted values of X12, yi2, and ^f^, with ler parabolic 
errors, are shown in the last column of Table II. In par- 
ticular, we obtain 



k f 2 [rad] = 0.07 ± 0.08 (la) 



(104) 



for parabolic errors. This is fully consistent (within la) 
with Eq. (|70[) for no new weak phases in decay, as ex- 
pected. To ascertain the impact on models in which new 
weak phases in decay are possible, we note that the re- 
lation between (f>® 2 and M^ p in Eq. (|T3"|) is modified, to 
very good approximation, as 



sin 0f 2 



MNP 



M 12 



sin 



-NP 



,SM 



)- 



(105) 



Therefore, we obtain the approximate (parabolic) la CL 
interval 



Im(M^) 
iMnl 



G [-0.01, +0.15] 



(106) 
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(for the usual phase convention in which Aff 2 M = is 
real), up to small corrections of a few percent or less from 
8(j)r/2, This is consistent, within lrr, with Eq. (|74| for 
no new weak phases allowed. Similarly, the correspond- 
ing HFAG (non-parabolic error) analysis, i.e., a direct fit 
to the experimental data which allows Ad ^ and in- 
corporates Eq. (|52|) . should be consistent with Eq. (|75|) . 

What will the sensitivity to new weak phases in de- 
cay be at a high luminosity flavor factory? We have 
seen that their impact on the time-dependent CP asym- 
metries (AY/) in SCS decays can be at most a few 
percent of the current errors (for D° — > n + ir~ and 
D° -> K+K-). In the case of CF and DCS decays (e.g., 
D° — > i^ s 7r + 7r~, i4T ± 7r T ) we saw that their current sensi- 
tivity to <fi is roughly an order of magnitude weaker than 
the maximal shifts allowed in (p^, and similarly for the 
precision with which Eq. (fBTj) . relating (y' + — y'~) and 
agL, can be tested. Thus, even with an order of magni- 
tude reduction in the errors on CPVINT and CPVMIX, 
as might be expected at a super- B factory with 75 ab _1 , 
it could be difficult to detect new weak phases in decay 
at currently allowed levels via time-dependent CP asym- 
metry measurements. 

The effects of new weak phases in decay are much easier 
to observe in D° and direct CP asymmetry measure- 
ments (in D° decays this requires comparison of the time- 
integrated and time-dependent CP asymmetries), as they 
are not suppressed by mixing (x or y). In particular, 
there is a good chance of detecting direct CP violation 
in SCS decays at a super-B factory, even if due solely to 
SM penguins. As a further illustration, we observe that 
the sum and difference of CP-conjugate time-integrated 
CP asymmetries in SCS decays would satisfy 



^L+R f a 7 



■^- + R f 4-a SL 



lfU f 

,<i 



y 



■ cos A 



R f aj 



a 



J- - Era?? 



R 



(107) 



up to negligible corrections of O(xr^j) and 0(yr 



' fj> 



Violations of the r/ = r j = relations satisfied by a/ and 
aj in Eqs. (|59|) and ([60]) could, therefore, be observed at 
a high luminosity flavor factory (e.g., in D° — > K*K de- 
cays) at currently allowed levels, unlike the violations of 
the corresponding (AY/ , AYj) time-dependent CP asym- 
metry relations in Eqs. ([92]) and (|93|) . 

Finally, it is interesting to observe that with sufficient 
statistics it could be possible to isolate and measure new 
contributions to arg(r 12 /r| 2 ), precisely because the ex- 
perimental sensitivity to (j) m individual decays substan- 
tially lags the direct CP asymmetry errors. Presumably, 
when the effects of non-zero values for (S(j>\ _ — Scpr/2) 
in Eqs. ([80 ]) -([84 ]) are observed, the direct CP asymmetry 
errors will be much smaller. If, for example, it turns out 
that 



V2|»K| 



(108) 



in the case of decay to a CP eigenstate, then we will 

have obtained a measurement of S4>p (since S<f>\ f ~ ai). 

For example, this situation could be realized: (i) in SCS 

decays, if new phases in CF/DCS amplitudes are near the 

current direct CP violation bounds, or (ii) in CF/DCS 

decays, in the more likely possibility that new phases 

only appear in SCS amplitudes. The determination of 

5(f>r could be combined with a measurement of (f>® 2 to fix 
lA/NPlsin^NP 



iff) in Eq. IT051) . 



C. B s — B s mixing 

Moving to B s — B s mixing, the CDF and DO collabo- 
rations are probing Sf, in B s — > and osl with 
combined uncertainties of 0.4 and 0.009, respectively 
Q. At LHCb with 2fb _1 the expected uncertainties are 
SS f w 0.02 H3 and <5a SL « 0.002 [||. If new subleading 
weak phases appear in B s decays to CP eigenstates, then 



Sf = r)f sign(a;) sin <j> Xf 



(109) 



with ipxj given in Eq. (|T6[) . The modified relation be- 
tween agL and Sf for decays to CP eigenstates [see Eq. 
([55]) ] follows from Eq. ([52]) . and is given to lowest order 
in rf and \q/p\ — 1 by 



7 1? P (1 _ ^2)1/2 =-sign(ycos<£ A/ ) 



sign(a; cos 4>\ f ) 
cos 2 <j)\ f 



Qsl 
2 



— ^— — 8(j>\. .(110) 



The new physics satisfies the inequality 
|(2r/ cos 5f sin <frf) tan <f>\ < 1 (unless <fi « ir/2), which 
implies that sign(y cos <f>\, ) = sign(y cos <f>) = +1 in the 
first term, and sign(x cos (f>\, ) = sign(a: cos <$) in the 
second term. 

There are no direct CP asymmetry measurements yet 
for B s decays mediated by b — > ccs transitions. How- 
ever, their magnitudes should be of same order as those 
for b — > ccs transitions in Bd decays. The best bound is 
ps 2%, from the direct CP asymmetry for B d — > J/^K° 
@ . As previously noted, the strong phase differences en- 
ter as sin Sf in the direct CP asymmetries, and as cos Sf 
in Stp\j and S(f>r- However, in the Standard Model the 
B — > J/^cf) amplitude is given by a particular color- 
suppressed combination of two effective operator ma- 
trix elements (Qi.2)- Moreover, the soft gluon contribu- 
tions to these matrix elements are for mally suppressed by 
AQCD/("ifcas) rather than Aqcd/"1{, [2^]. We also note 
that significant strong phase differences (~ 30° — 50°) 
between the different isospin amplitudes in B — ► D^n 
and B — ► D^*'K decays are known to exist, due to color- 
suppressed channels [27]. Thus, significant strong phase 
differences between the SM and new physics B s —> J/fycj) 
or Bd — > J/^K amplitudes can be expected. We there- 
fore take 



*A,|, l^rl <5%, 



(111) 
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for new physics in b — » ccs transitions. Effects of this size 
in the second term on the r.h.s. of Eq. (|110[) . applied to 
B s — > J would be difficult to observe at LHCb, given 
an order of magnitude larger projected experimental un- 
certainty on the first term 



0(0.4) 



(112) 



for 2 fb 1 (obtained from Sash above, and the SM central 
value for \x/y\ in [H). 

New CP violating effects at the 5% level in the tree- 
amplitudes would be quite exotic. If the new physics 
enters the b — > ccs transitions via gluonic or electroweak 
penguins, which we believe is a far more likely scenario, 
then its contributions to S<fi\ J/<!!<t> and S(fir would be neg- 
ligible. Recall that new CP violating amplitudes in pen- 
guin dominated Bd decays, e.g., Bd — > (fiK s , are con- 
strained to lie below 0(10%). 

Finally, 5(fir receives a significant Standard Model con- 
tribution (relative to the leading [A c /A*] 2 CKM structure 
in T 12 /T* 12 ) 



SM. 



— 41m 



1 12 

F 



12 / SM 



1 12 
pec 
1 12 



(113) 



SM 



where Xi = V* s Vib- and are defined in the Stan- 
dard Model expression for Ti2, 

1 12 — -\ A c 1 1 



9 \ \ p«c 



+ \ 2 u r™). (114) 



In the OPE treatment they differ only with respect to 
quark content in loops: two charm quarks vs. one charm 
and one up quark, and satisfy — T"| (28[. This is 
likely to be the dominant contribution to Stfir, certainly 
if new physics only enters through gluonic or electroweak 
penguins. With sufficient statistics it could be possi- 
ble to isolate and measure 6 (fir via Eq. (|110[) applied 
to B s — » J/^fcfi. This would require that the hierarchy in 
Eq. (|108[) , equivalent to |<5</>r/2| ~> \S(fi\ J/qjit> \, is satisfied. 
In practice, a substantial improvement of the direct CP 
asymmetry bounds for b — * ccs transitions would also be 
required. 



VII. DISCUSSION AND CONCLUSION 

If 012 = arg(Mi2/ r i 2 ) is the only source of CP vi- 
olation in neutral meson decays, then CP violation in 
pure mixing (CPVMIX), i.e., \q/p\ ^ 1, and CP viola- 
tion in the interference of decays with and without mix- 
ing (CPVINT), i.e., <fi 7^ 0, are related phenomena. More- 
over, (fi would be related to the underlying mixing param- 
eters | A'/i2 1 , | Ti2 1 , and <fi\2 of relevance to model building. 
New weak phases in the decay amplitudes would enter 
and modify these relations. However, existing direct CP 
violation measurements provide stringent constraints on 
their magnitudes in the (tree-level dominated) D° and 
B s decays of interest to us, implying that any modifica- 
tions to the relations between CPVMIX and CPVINT 



must be small perturbations. We summarize these re- 
sults, and their implications below. 

The general relation between <fi and \q/p\ (CPVINT 
and CPVMIX) in the limit of no new weak phases in de- 
cay is derived in Section IV, see Eq. (|54|) . It leads to cor- 
relations between the semileptonic and time-dependent 
CP asymmetries and additionally, in the D° system, be- 
tween the semileptonic and time-integrated CP asymme- 
tries. We remind the reader that in D a decays the time- 
dependent (AY/) and time-integrated (a/) CP asymme- 
tries must be equal in the limit of no direct CP violation 
[ill ] (no new weak phases in decay), see Section III. 

Below we will refer to the whole complex of re- 
lations obtained via applications of Eq. (|5"4")) as the 
CPVMIX/CPVINT relations. We give them a fairly 
complete treatment in the case of tree-level dominated 
D° decays, covering singly Cabibbo suppressed (SCS) 
decays to CP {K+R-, tt+tt') and non-CP (K*K) 
eigenstates, as well as Cabibbo favored (CF) and dou- 
bly Cabibbo suppressed (DCS) decays to CP eigen- 
states (K s ir°) and to "wrong-sign" non-CP eigenstates 
(K + tt~), where examples are included in parentheses. 
In the case of B s decays, we confirm the correlation be- 
tween the semileptonic and time-dependent CP asymme- 
tries obtained in 0] for decays to CP eigenstates, and we 
also obtain the correlation for decays to non-CP eigen- 
states, see Eqs. (frJ5)) and 

For SCS D° decays to non-CP eigenstates, CP con- 
jugate decay rates are of same order. Therefore, pairs 
of time-dependent and time-integrated CP asymmetries 
are accessible to experiment. We find that the relation 
between <fi and \q/p\ implies that the ratio of CP asymme- 
tries within each pair is given by the inverse ratio of the 
corresponding CP averaged decay rates, see Eq. By 
contrast, the near equality of CP conjugate pairs of time- 
dependent CP asymmetries (Sf and Sj) for B s decays to 
non-CP eigenstates is a trivial consequence of y -C x and 

The general expression derived for (fi in terms of the 
underlying mixing parameters I-M12I, and <fi\2 1 in 

the limit of no new weak phases in decay, is given in Eq. 
([52")) . It can be combined with similar expressions for 
x, y, and \q/p\, see Eqs. (fl"4l) . (fT5|) . and ([4"8)l. to extract 
the underlying D a — D° mixing parameters |M^|, |rf^|, 
and (fii2 from a direct fit to the experimental data. In 
this work we adopt the simpler strategy of extracting 
the mixing parameters from a fit to the HFAG outputs, 
which include x, y, \q/p\, and (fi. The (parabolic) HFAG 
output error matrix is used to construct a \ 2 function, 
see Section V. We find that 

• (fii2 is currently being probed at the level of 0.10 
[rad] at ler, see Eq. ([701) . 

• Incorporating the relation between cfi and (^12, Eq. 
(|52p , into the fit reduces the experimental errors on 
(fii 2 by a factor of three for the current data set. 

The preliminary (non-parabolic) HFAG fit directly to the 
data [3, also obtained using Eqs. 1(131). (USD, l[35]). and 
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(|52f for x, y, \q/p\, and (j>, yields a sensitivity to <f>^ 2 °f 
0.10 [rad] (la); the 95% CL bound is 0f 2 < 0.30 [rad], 
see Eq. (|7Tjl . 

Two questions concerning the impact of new weak 
phases in decay need to be addressed: (i) to what ex- 
tent can the CPVMIX/CPVINT relations be violated in 
the D° and B s systems, and how well could such viola- 
tions be measured in the future; (ii) to what extent can 
the bound on <p^ 2 be modified. The violations in (i) can 
be characterized precisely in terms of the shifts 54>\, in 
the CPVINT observables 4>\ f = &rg\qAf /pAf\ with re- 
spect to 4>, the shift 6<j>r in argfT^/r^) with respect to 
the appropriate leading SM contribution, and the direct 
CP asymmetries a d f (see Eqs. (|89]) - ([93]) . Eqs. lp |l -([T00 l) . 
and Eq. Q110p in Section VI) . Thus, we need to know how 
large these quantities can be. 

Direct CP violation bounds provide stringent con- 
straints on subleading amplitudes containing new weak 
phases. Strong phase differences enter as cos S in the 
direct CP asymmetries, and as sin 5 in 6<f)\ f and 54>v- 
However, we argue that in all cases of interest the new 
physics amplitudes would have significant strong phase 
differences with respect to the leading SM amplitudes 
(due, essentially, to different color and chirality struc- 
tures for the underlying effective operators). Therefore, 
the direct and time-integrated CP violation bounds also 
translate into order of magnitude bounds on 5<fi\,, S(p\^, 
and 54>r due to new weak phases in decay (see Eqs. (|88[). 
([97]) . and ifTTT]) for the SCS, CF/DCS, and b -» ccs tran- 
sitions, respectively). 

The main implications of these bounds for D° — D a 
and B s — B s mixing today are: 

• In SCS D° decays the maximal allowed viola- 
tions of the CPVMIX/CPVINT relations are of 
0(a few%) of the current experimental errors on 
the time-dependent CP asymmetries (CPVINT), 
see Eqs. (|59 ]) -([9"2 |) . 

• In CF/DCS D° decays the maximal allowed vio- 
lations of the CPVMIX/CPVINT relations are of 
0(10%) of the current experimental errors on the 
time-dependent CP asymmetries (CPVINT), see 

Eqs. pHing]). 

• Violations of Eq. (|52|) . relating <f> and 4>® 2 , are sim- 
ilarly bounded relative to the present day exper- 
imental sensitivity to in SCS and CF/DCS D° 
decay modes, respectively. 

• Consequently, the bounds on <f>^ 2 m Eqs. (|7D|) and 
(I7ip can not be significantly modified by new weak 
phases in decay, see Eq. (|104[) . 

• For b — > ccs transitions, the maximal allowed vi- 
olation of the B s CPVMIX/CPVINT relations is 
0(5%) in absolute terms, see Eq. (jllOp . 

At a high luminosity flavor factory (with 75 ab -1 ), we 
assume that there will be an order of magnitude improve- 



ment in precision for individual D° — D° mixing measure- 
ments, and in the global fit to the data (a reduction of 
sa 6 in the errors on <j> and \q/p\ for the global HFAG fit 
is projected in [29]). The error on the B s semileptonic 
CP asymmetry at LHCb (with 2fb~ 1 ) is expected to be 
(5osl ~ 0.002 [23. Therefore, our conclusions on the sen- 
sitivity of mixing measurements at these facilities to new 
weak phases in decay are: 

• Violations of the D° - W CPVMIX/CPVINT re- 
lations will be probed at the same order as the cur- 
rently allowed maximal violations (obtained from 
direct and time-integrated CP violation measure- 
ments), implying that they could be difficult to ob- 
serve at a high luminosity flavor factory. 

• The "goodness" of a global fit to the D° -W mix- 
ing data which assumes no new weak phases in 
decay would probably be more sensitive to their 
effects than violations of the CPVMIX/CPVINT 
relations in individual decay modes. 

• The expected error in the B s semileptonic CP 
asymmetry at LHCb is prohibitively large for a 
meaningful probe of the B s CPVMIX/CPVINT re- 
lations to be carried out. 

In principle, with sufficient statistics it would be pos- 
sible to determine 64>r in the D° and B s systems, if the 
violations of the CPVMIX/CPVINT relations are much 
larger than the direct CP asymmetries in the SCS or 
CF/DCS transitions and the b — > ccs transitions, respec- 
tively (see the discussions at the ends of Sections VIB 
and VIC). 

We emphasize that the D° and D direct CP asym- 
metry measurements provide much more sensitive probes 
of new weak phases in decay than the time-dependent 
CP asymmetries (which correspond to differences be- 
tween the D° and D° time of decay profiles). Recall 
that in D° decays, the direct CP asymmetries are ob- 
tained from comparison of the time-integrated and time- 
dependent CP asymmetries. The effects of new weak 
phases in the time-dependent CP asymmetries are nec- 
essarily suppressed by mixing (x or y). Therefore, the 
most likely scenario at a high luminosity flavor factory is 
that improved precision in the time-integrated or direct 
CP violation measurements will imply that the effects of 
new weak phases in decay lie beyond the reach of the 
time-dependent CP asymmetry measurements. It will of 
course still be possible to probe for new weak phases in 
b — ► ccs transitions at a super-B factory (B decays) and 
at the LHC (B and B s decays) via direct CP violation 
measurements . 

Finally, and of immediate interest for CP violation 
in D° — D° mixing, the bounds on <p^ 2 imply that 
Im(M^ p )/|Mi 2 | is being probed at the 10% level at la 
(for the usual phase convention in which Mf^ is real, 
and where M^ 2 P is the new physics contribution). The 
preliminary HFAG 95% CL interval for </>f^ implies that 
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\lm(M^)/M 12 \ < 0.30 at 95% CL. These results apply 
to models without new weak phases in decay, or with 
new weak phases in decay (up to an additive correction 
of less than a few percent), see Eqs. ([73" )) ([75]) . and Eqs. 
1(105)1 . (fTOo) . Examples in which CP violation in D° - IP 
mixing at such levels is possible have recently been dis- 
cussed in the context of supersymmetry, little Higgs mod- 
els, warped extra dimension models, and the minimal fla- 
vor violation framework @, [ljl [3(| [HI, [HJ . 



Acknowledgements. We thank Yuval Grossman and 
especially Yossi Nir for numerous discussions on the mix- 
ing and CP violation formalism, Gilad Perez for a discus- 
sion of the bound on CP violation in D° — D° mixing, 
and Alan Schwartz for providing us with the HFAG out- 
put correlations (for Ad — 0), the preliminary HFAG fit 
results for 4>i2, and discussions. A. L. K. is supported by 
DOE grant FG02-84-ER40153, M. D. S. by NSF grant 
PHY0757876. 



[1] Y. Nir, SLAC-PUB-5874, Lectures given at 20th Annual 
SLAC Summer Institute on Particle Physics (1992). 
Y. Nir, |arX~iv:hep-ph/0510413| 

G. Valencia, Phys. Rev. D 39, 3339 (1989); A. Datta 
and D. London, Int. J. Mod. Phys. A 19, 2505 (2004) 
|arXiv:hep-ph70 303159 . 

S. Bergmann, Y. Grossman, Z. Ligeti, Y. Nir and 
A. A. Petrov, Ph ys. Lett. B 486, 418 (2000) 

|arXiv:hep-ph70 005181 . 

Y. Grossman, Y. Nir and G. Perez. larXTv:0904.0305 1 [hep- 
ph]. 

I. I. Bigi, M. Blanke, A. J. Buras and S. Recksiegel, 
larXiv:0904.1545l [hep-ph] . 

M. Ciuchini, E. Franco, D. Guadagnoli, V. Lubicz, 
M. Pierini, V. Porretti and L. Silvestrini, Phys. Lett. B 
655, 162 (2007) [arXiv:hep-ph /0703204| . 
J. M. Link et al. , [FOCUS Collaboratio n], Phys. Lett. B 
485, 62 (2000) [arXiv:hep-ex/0004034] . 
E. Barberio el al, [HFAG Collabora- 
tion], arXiv:0808.1297 and updates at 



ht t p : / / www . slac . Stanford . edu/xorg/hfag 

C. Amsleret al., [PPG C ollaboration] , Phys. L ett. B 667, 

1 (2008) and updates at |http://pdg.lbl.gov/| 

Y. Grossman, A. L. Kagan and Y. Nir, Phys. Rev. D 75, 

036008 (2007) [arXiv:hep-ph/0609178| . 

J. L. R osner and D. A. Suprun , Phys. Rev. D 68, 054010 

(2003) [arXiv:hep-ph/0303117l . 

L. M. Zhang et al. [BE LLE "Col laboration] , Phys. Rev. 
Lett. 96, 151801 (2006) [arXiv:hep-ex/0601029l . 
Alan Schwartz, private communication 
O. Gedalia, Y. Grossman, Y. Nir and G. Perez, 
larXiv:0906.1879l [hep-ph] . 

A. J. Schwartz, talk at Charm 2009, Leimen, May 2009. 
M. Bobrowski, A. Lenz, J. Riedl and J. Rohrwild, 



|arXiv:0904.3971 [hep-ph]. 
[18] M. Staric et al. [B elle Co llaborati on], Ph ys. Rev. Lett. 

98, 211803 (2007) [arXivihep ex/0703036| ; B. Aubert et 
al. [BABAR Collaboration], Phys. Rev. D 78, 011105 
(2008) larXiv:0712. 2249 [hep-ex]]. 

[19] A. L. Kagan, arXiv:hep-ph/0407076 

[20] B. Aubert et al. [BABA R Collaboration], Ph ys. Rev. 

Lett. 98, 211802 (2007) [arXiv:hep-ex/0703020] . 
[21] S. Bergma nn and Y Nir, JHEP 9909, 031 (1999) 

[arXiv:hep-ph/990939l] . 
[22] G. D'A mbrosio and D. N. Gao , Phys. Lett. B 513, 123 

(2001) [arXiv:hep-ph/0105078| . 
[23] K. Abe et al. [BELLE Collaboration], Phys. Rev. Lett. 

99, 131803 (2007) |arXiv:0704.1000l [hep-ex]]. 

[24] M. N. Minard, talk at "The Year of the Ox", Aspen Win- 
ter Conference, February 2009. 
[25] N. Brook et al., LHCb-note 2007-054. 
[26] M. Beneke, G. Buchalla, M. Neubert and 
O T. Sachrajda, Nu cl. Phys. B 591, 313 (2000) 
[arXiv:hep-ph/0006124] . 

S. Mantry, D. Pirjol and I. W. Stewart , Phys. Rev. D 68, 
114009 (2003) [arXiv:hep-ph/0306254] . 
A. Lenz and U. Nier ste, JHEP 0706, 072 (2007) 



[27 
[28 
[29 
[30 
[31 
[32 



[arXiv:hep-ph/0612167 . 

A. J. Schwartz, plenary talk at CIPANP 2009, San Diego, 
May 2009. 



K. Blum, Y. Grossman 
larXiv:0903.2"Tl"8l [hep-ph] . 
C Csaki, G. Perez, Z. 
larXiv:0907.0474l [hep-ph] . 
A. L. Kagan, G. Perez, 
larXiv:0903.i794l [hep-ph]. 



Y. Nir and G. Perez, 
. Surujon and A. Weiler, 
T. Volansky and J. Zupan, 



